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ABSTRACT 



The phase structures of Nambu-Jona-Lasinio models with one or two flavours have been in- 
vestigated at non-zero values of fi and H, where H is an external magnetic field and fi is the 
chemical potential. In the phase portraits of both models there arise infinitely many massless 
chirally symmetric phases, as well as massive ones with spontaneously broken chiral invariance, 
reflecting the existence of infinitely many Landau levels. Phase transitions of first and sec- 
ond orders and a lot of tricritical points have been shown to exist in phase diagrams. In the 
massless case, such a phase structure leads unavoidably to the standard van Alphen-de Haas 
magnetic oscillations of some thermodynamical quantities, including magnetization, pressure 
and particle density. In the massive case we have found an oscillating behaviour not only for 
thermodynamical quantities, but also for a dynamical quantity as the quark mass. Besides, in 
this case we have non-standard, i.e. non-periodic, magnetic oscillations, since the frequency of 
oscillations is an if-dependent quantity. 
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1. Introduction 



The exploration of strongly interacting matter at high density and in the presence of external 
electromagnetic fields is of fundamental interest and has potential applications to the quark- 
gluon plasma and heavy-ion collisions, to cosmology and astrophysics of neutron stars. Recently, 
some aspects of this problem were considered in where it was pointed out that in QCD at 
high density a new phase with colour superconductivity might exist. The influence of external 
magnetic fields on the QCD vacuum was, for example, studied in 0. 

Our goal is to investigate properties of the strongly interacting cold quark matter in the 
presence of both external magnetic field H and non-zero chemical potential /i. The subject 
is closely related to magnetic oscillations of different physical quantities. In this connection 
we should remember that the van Alphen-de Haas effect (oscillations of the magnetization) 
was for the first time predicted by Landau and then experimentally observed in some non- 
relativistic systems (in metals) more than sixty years ago 0, ^. At present, a lot of the 
attention of researchers dealing with magnetic oscillations is focused on relativistic condensed 
matter systems (mainly on QED at /x, if 7^ 0), since the results of these studies may be applied 
to cosmology, astrophysics and high-energy physics P, ^ . 

It is well known that up to now the consideration of QCD at /i, if 7^ is a difficult problem. 
This is partly due to the fact that numerical lattice simulations at /i 7^ have not been able to 
overcome problems associated with the complex part of the fermionic determinant. Moreover, 
the incorporation of a magnetic field into lattice gauge calculations is not elaborated sufficiently, 
either. For these reasons, when considering quark matter at /i, if 7^ 0, many authors prefer to 
deal with adequate models (e.g. with the MIT bag model [0), rather than with QCD. 

In the present paper we shall study the above problem in the framework of some specific 
QCD-like quark models. Namely, we shall investigate the influence of an external magnetic 
field and chemical potential on the vacuum structure of Nambu-Jona-Lasinio (NJL) models 
containing four-fermion interactions The simplest one, denoted as model I, refers to the 

one-flavour case and is presented by the Lagrangian 

G 

Li = Qkidqt + ^[{mkY + {qkiimk)\ (1) 

where all quark fields belong to the fundamental multiplet of the colour SU{Nc) group (here 
the summation over the colour index k = 1, ...,Nc is implied). Obviously, Li is invariant under 
(global) SU{Nc) and U{l)v transformations as well as continuous U{1)a chiral transformations 

Qk^e^'^^qk ; {k=l,...,N,). (2) 

The second, more realistic case considered here and referred to as model II is a two-flavour 
NJL model whose Lagrangian has the form 

G 

L2 = qidq+-—[{qqf + {qi^''fq)\ (3) 
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where g is a flavour isodoublet and colour- A^c-plet quark field and r are isospin Pauli matrices 
(in (3) and below, flavour and colour indices of the quark field q are now suppressed). The 
Lagrangian L2 is invariant under (global) U{2) f x SU{Nc) as well as under chiral U{2)iX U{2)r 
groups. 

NJL models were proposed as a good laboratory for investigating the non-perturbative 
phenomenon of dynamical chiral symmetry breaking (DCSB), which occurs in the physics of 
strong interactions, as well as for describing the low-energy sector of QCD (see e.g. papers 
[|lO| , [TI| , [1^ and references therein). Since there are no closed physical systems in nature, the 
influence of different external factors on the DCSB mechanism is of great interest. In this 
relation, special attention has been paid to the analysis of the vacuum structure of NJL-type 
models at non-zero temperature and chemical potential in the presence of external 

(chromo-)magnetic fields |]15|, |T6|, 0, with allowance for curvature and non-trivial space-time 
topology |TS|, |TP|. The combined influence of external electromagnetic and gravitational fields 
on the DCSB effect in four-fermion field theories was investigated in 
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In the present paper the phase structures and related oscillating effects of the above- 
mentioned Nambu- Jona-Lasinio models are considered at ii,H ^ 0. We will show that, here, 
the set of oscillating physical parameters in NJL models is richer than in QED at /i, if 7^ 0. Be- 
sides, in the NJL models, in contrast to QED and similar to some condensed- matter materials, 
there exist non-periodic magnetic oscillations. 



2. NJL models at ^ and H = 

First of all let us prepare the basis for the investigations in the following sections and consider 
in detail the phase structure of the model I at non-zero chemical potential 7^ and H = 0. 

Recall some well-known vacuum properties of the theory (1) at /i = 0. The introduction of 
an intermediate quark-meson Lagrangian 

Li = qidq - g((Ti + i(T2l5)q ~ ^(^1 + ^2) (4) 

greatly facilitates the problem under consideration. (In (4) and other formulae below we have 
omitted the fermionic index k for simplicity.) Clearly, using the equations of motion for the 
bosonic fields ai^2, the theory in (4) is equivalent to that in (1). From (4) we obtain the one-loop 
expression for the effective action 



exp^iNcSef f{(Ti^2)) = J DqDqexp J Lid'^x 



where 



5'e//(o-i,2) = - j ^'^^ ^^2G^^ -i\ndet{id - ai - nso-a). 
Assuming that in this formula ai^2 are independent of space-time points, we have by definition 

'S'e//(0-i,2) = -Vo(cri,2) J d'^X, 
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Vo{<yi,2) = § + 2^ / - P") = ^o(S), (5) 

where S = \J crl + 02. Next, introducing in (5) Euclidean metric (po — ^ Wo) and cutting off the 
range of integration {p^ < A^), we obtain 

The stationarity equation for the effective potential (6) has the form 

Now one can easily see that at G < Gc = 47r^/A^ Eq. (7) has no solutions apart from S = 0. 
Hence, in this case fermions are massless, and chiral invariance (2) is not broken. 

If G > Gc, then Eq. (7) has one non-trivial solution, So(G, A) 7^ 0, such that F(Eo) = 0. In 
this case Sq is a point of global minimum for the potential Vo(E). This means that spontaneous 
breaking of the symmetry (2) takes place. Moreover, fermions acquire a mass M = So(G, A). 

Let us now consider the case where /i > and the temperature T 7^ 0. In this case, the 
effective potential V^t{^) can be found if the measure of integration in (5) is transformed in 
the standard way according to the rule p2 



Summing there over n and letting the temperature in the obtained expression tend to zero, we 
obtain 

V,{T.) = Vo(S) - 2 J - - ^S^), 

where 6{x) is the step function. Finally, by performing the momentum integration, we find 



2/iV/i2-S2 + S^ In 



It follows from (8) that, in the case G < Gc and at arbitrary values of the chemical potential, the 
chiral symmetry (2) is not broken. However, at G > Gc the model has a rich phase structure, 
which is presented in Fig. 1 in terms of fi and M. (At G > Gc one can use the fermionic mass 
M as an independent parameter of the theory. The three quantities G, M and A are connected 
by Eq. (7).) In this figure the solid and dashed lines represent the critical curves of the second- 
and first-order phase transitions, respectively. Furthermore, there are two tricritical points 0, 

point of the phase diagram is called a tricritical one if, in an arbitrarily small vicinity of it, there are 
first- as well as second-order phase transitions. 
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a and /3, two massive phases B and C with spontaneously broken chiral invariance as well as 
the symmetric massless phase A in the phase portrait of the NJL model I (detailed calculations 
of the vacuum structure of this NJL model can be found in |]14|)- In the present model the 



dynamical quark mass So(/i), given by the global minimum of the potential V^(S) as a function 
of /i behaves as depicted in Fig. 2. 

We should note that the phase transition from i? to C is the quark-matter analogue of the 
so-called insulator-metal phase transition in condensed-matter physics. This is due to the fact 
that, in the vacuum of phase the particle density (the analogue of conductivity electron 
density in condensed materials) is zero, while in the vacuum of phase C there arises a non-zero 
density of charged particles, so that it looks like a Fermi-liquid ground state of metals. 

To investigate the vacuum properties of the two-flavour NJL model (model II) it is again 
convenient to employ, instead of the quark Lagrangian (|]), the equivalent quark-meson La- 
grangian 

U = qtYd^q - q{cr + ^7'r7r)g - ^(^' + ^')- (9) 

Using here calculations similar to the case with model I, one can find the effective potentials 
for /i = and /x 7^ expressed in terms of meson fields a, vf . These potentials have the form 
of the potentials (6) and (8) for the model I, respectively, with the exception that the factor 
(IGtt^)^^ in (6), (8) has to be replaced by (87r^)~^. Moreover 0, in the case under consideration 
T? = cr^ + TT^. It follows from this similarity that the phase structure and the phase portrait of 
the model II are qualitatively the same as those of the model I (see Fig. 1). 



3. Phase structure of the model I at /i 7^ and H ^ 

In the present section we shall study vacuum magnetic properties of NJL systems. For the 
model I at /i = 0, this problem was considered in |jl5|, [l^. It was shown in [jl5l that at G > Gc 
the chiral symmetry is spontaneously broken for arbitrary values of the external magnetic field 
H, including the case H = 0. At G < Gc the NJL model I has a symmetric vacuum at H = 0. 
However, if an external (arbitrarily small) magnetic field is switched on, then for all G G (0, Gc) 
there is a spontaneous breaking of the initial U{1)a symmetry (2) |jT^. This is the so-called effect 
of dynamical chiral-symmetry breaking (DCSB) catalysis by an external magnetic field. (This 
effect was observed for the first time in the framework of a (2+l)-dimensional Gross-Neveu 
model in and was then explained in Now this effect is under intensive investigations 



since it has a wide range of possible applications in physics^.) 

Let us recall some aspects of the problem at /i = 0, 7^ 0. In order to find in this case the 
effective potential Vff(S) of the NJL model I, gauged by an external magnetic field according 



to = — ieA^, = 5^2XiH, one can use the well-known proper-time method |26 



^Chiral symmetry of the two-flavour NJL model (3) is realized in the space of meson fields (cr, tt) as the 
rotation group 0(4), which leaves the "length" of the particle vector (ct, tt) invariant. Hence, all effective 
potentials of this model depend on the single variable E = \/a^ + 7?^. In the following wc choose tt = 0, 
assuming the absence of a pion condensate. 

^Some recent references on this subject are presented in psj. 
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or momentum-space calculations [27|, which give the following expression: 

MS) = + g + II f |exp(-,E Voth(e/f.). 

In this formula e has a positive value. It is useful to rearrange this expression in the form 

Vh(S) = + + (10) 



where 



H2 (g^)2 



(eiJs) coth(eiJs) - 1 - i^E^ 



(11) 



thereby isolating the contributions of the matter, the field and the electromagnetic interaction 
energy densities exphcitly. The potential Vo(S) in (|1^) is up to an unimportant (infinite) 
additive constant not depending on E, equal to expression (5). Hence, the ultra-violet (UV)- 
regularized expression for it looks like (6). 

The integral of the function Z^E) is also UV- diver gent, so we need to regularize it. The 
simplest possibility is to cut it off at the lower boundary, which yields 



7 being the Euler constant. Clearly, the last term in (|T2D contributes to the renormalization of 
the magnetic field and electric charge, in a way similar to what occurs in quantum electrody- 



namics 



The potential Vff(S) in (|TT]) has no UV divergences, so it is easily calculated with the help 
of integral tables p8[|. The final expression for Vff(S) in terms of renormalized quantities is 
then given by 



Vh{^) = ^ + Vo{E) - ^{C'(-l, x) -\[x'-x]\nx+ ^}, (13) 



where x = T?/{2eH), C^{v,x) is the generalized Riemann zeta-function and C,'{—l.,x) = 
dC,{v.,x) / dv\v=-i. The global minimum point of this function is the solution of the stationarity 
equation 

^^h(S) = ^{F(S)-/(S)} = 0, (14) 
where F{Y?) is given in (7), and 

/(S) = 2e//{lnr(x) - ^ln(27r) +x- ^(2x- l)lnx} 

= / — exp(-sS^)[ ei/scoth(ei/s) - 1]. (15) 
Jo 



One can easily see that there exists, for arbitrary fixed values of H, G, only one non-trivial 
solution I1q{H) of Eq. (0), which is the global minimum point of Vh{^)- There, Tiq{H) is a 
monotonically increasing function of H, and at — *• cxd 



2o(ff) 



vr V 12 



(16) 



However, at H 



eH 



A^)}, ifG<a, 
if G > Gr. 



(17) 



So, at G < Gc and H = the NJL vacuum is chirally symmetric, but an arbitrarily small 
value of the external magnetic field H induces DCSB, and fermions acquire a non-zero mass 
So(-ff) (the magnetic catalysis effect of DCSB). 

Now let us consider the more general case, when H ^ and fi ^ 0. In one of our previous 
papers an effective potential of a 3D Gross-Neveu model at non-zero H, /i and T was 
obtained. Similarly, one can find the effective potential in the NJL model I at H,T, fi ^ 0: 



V, 



TeH 



k=0 



exp 



-/3(£fc+Ai) 



1 + exp" 



-/3(£fc-At) 



(18) 



where /? = l/T, = 2 — 6ok, £k = V^'^ + p'^ + 2eHk, with k = 0, 1,2, ... denoting Landau 
levels, and the function Vff(E) is given in (|13|). If we let the temperature in (p!8D tend to zero, 
we obtain the effective potential of the NJL model I at if, /i 7^ 0: 



47r J -00 



(19) 



which, by performing the integration, can easily be cast into the form 



V, 



Hfj.\ 



eH 

47^2 



E «fc^(^ - Sk) \ f^y fJ'^ - sl - si In 



k=0 



(20) 



where = -|- 2eHk. Finally, let us present the stationarity equation for the potential 



47r2 



F(E) - /(S) + 2eH ^ - Sk) In 

fc=0 



Sk 



0. 



(21) 



In order to get a phase portrait of the model under consideration we should find a one-to- 
one correspondence between points of the [fi, if )-plane and the global minimum points of the 



function (|20|) , i.e. by solving Eq. (|2l|) we should find the global minimum S(/i,if) of the 
potential (pO|) and then study its properties as a function of (/i, H). 
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3.1 The case G < G^- Magnetic catalysis and chemical potential 

In order to greatly simplify this problem, let us divide the (/i, i7)-plane into a set of regions Uk- 

oo 

= y cjfe, uJk = {{ii,H):2eHk<ii^<2eH{k+l)}. (22) 

k=0 

In the ujo region only the first term from the series in ( pO] ) and ( pT] ) is non-vanishing. So, one 
can find that, for the points (/i, H) G ujq, which are above the line ^={(yU, H) : fi = So(-ff)}, the 
global minimum is at the point S = 0. Just under the curve / the point S = So(i^) is a local 
minimum of the potential (|20|), whereas S = So(if) becomes a global minimum, when {n,H) 
lies under the curve /i = fic{H), which is defined by the following equation: 

Vh,{0) = Vh,{^o{H)). (23) 

Evidently, the line /i = fJ^dH) is the critical curve of first-order phase transitions. In the Uq 
region Eq. ( p3D is easily solved: 

97r 

/ic(i^) = -r^lVniO) - Vni^omf'. (24) 

Using the asymptotics ([T7|) of the solution T,o{H) at — > 0, we find the following behaviour 
of iXc{H) at if ^ 0: 

Hence, we have shown that at /i > fJ^dH) {G < Gc) there exists a massless symmetric phase 
of the NJL model (numerical investigations of ( pO]) and (plD give us a zero global minimum 
point for the potential Vff^(S) in other regions uji,uj2, ■■■ as well). The external magnetic field 
ceases to induce the DCSB at /i > fic{H) (or at sufficiently small values of the magnetic field 
H < Hc{fi), where H^fJ,) is the inverse function of Hc{H)). However, under the critical curve 
(|2^ (or at H > Hc{fi)), owing to the presence of an external magnetic field, the chiral symmetry 
is spontaneously broken. Here the magnetic field induces a dynamical fermion mass Ilo{H), 
which has a /i- independent value. 

Lastly, we should also remark that in the NJL model (1) the magnetic catalysis effect takes 
place only in the phase with zero particle density, i.e. at /i < fic{H). If /x > fic{H), we have 
the symmetric phase with non-zero particle density, but here the magnetic field cannot induce 
DCSB. 

3.2 The case G < Gc- Infinite cascade of massless phases 

In the previous subsection we have shown that the points (/x, H), lying above the critical curve 
/i = Hc{H), correspond to the chirally symmetric ground state of the NJL model. Fermionic 
excitations of this vacuum have zero masses. At first sight, it might seem that the properties of 
this symmetric vacuum are slightly varied, when parameters fi and H are changed. However, 
this is not the case, and in the region fi > ^c{H) we have infinitely many massless symmetric 
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phases of the theory corresponding to infinitely many Landau levels, as well as a variety of 
critical curves of second-order phase transitions. Let us next prove this. 

It is well known that the state of thermodynamic equilibrium (the ground state) of an 
arbitrary quantum system is described by the thermodynamic potential (TDP) Q, which is just 
the value of the effective potential at its global minimum point. In the case under consideration, 
the TDP H) at > fic{H) has the form 



eH 



47]-2 



VniO) - 

oo 



k=0 



In 



4 + 



)/«]} 



(25) 



where = yleHk. We shall use the following criterion of phase transitions: if at least one 
first (second) partial derivative of f2(/i, H) is a discontinuous function at some point, then this 
is a point of a first- (second-) order phase transition. 

Using this criterion, let us show that boundaries of Uk regions (^2]), i.e. lines Ik = {(yU, H) : 
fi = y/2eHk} [k = 1,2,...), are critical lines of second-order phase transitions. In an arbitrary 
Uk region the TDP (ESf) has the form: 



= Vh{0) - 
eH 



An- 




- ef In 



(26) 



From ( P6|) one easily finds 



as well as: 



on. 



d'^k-i 



0, 



2vrV/i^-4 



-oo. 



(27) 



(28) 



Equation (|27|) means that the first derivative dQ/dfi is a continuous function on all lines Ik- 
However, the second derivative d'^Q/{dn)'^ has an infinite jump on each line Ik (see (pS])), so 
these lines are critical curves of second-order phase transitions. (Similarly, we can prove the 
discontinuity of d'^Vt/^dHY' and d'^Vt/ djjidH on all lines In-) 

The results of the above investigations are presented in Fig. 3, where the phase portrait of 
the NJL model I at G < Gc in the (/i, if)-plane is displayed. 



3.3 The case G > Gc 

Concerning supercritical values of the coupling constant, we shall consider here only the case 
Gc < G < (L225...)Gc, where the phase portrait of the model I is quahtatively represented 
in Fig. 4. In this figure one can see infinite sets of symmetric massless Aq^Ai,--- phases, as 
well as massive phases Cq, Ci,... with DCSB. In addition, there is another massive phase B- 
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Dashed and solid lines in Fig. 4 are critical curves of first- and second-order phase transitions, 
respectively. One can also see on this portrait infinitely many tricritical points tk,Sk {k = 
0,1,2,...). For a fixed value of k the point tk lies inside, but the point Sk is on the left 
boundary of the corresponding Uk region p^). Each critical line Ik coincides with a part of the 
ujk boundary. In Table 1, we give the values of the external magnetic field corresponding to 
tricritical points to and sq. 

The presence of an infinite cascade of massless A^-phases in the case G > Gc may be proved 
in a way similar to what was done in the previous subsection. However, now an infinite set 
of massive chirally non-symmetric phases is available thanks to the particular structure of the 
function in (^). A detailed numerical investigation of this function shows that, for some 
values of (/x, H) inside the Ck region (see Fig. 4), as a function of S qualitatively behaves 
like the curve, drawn in the Fig. 5. At these values of {fi, H) there is only one non-trivial 
solution HkilJ-jH) of the stationary equation PT|), which is the global minimum point of the 
effective potential (|20D and at the same time is the quark mass in the phase Gk of the theory. 
Remark that in each phase Gk the quark mass is a //-dependent function. In contrast, in the 
phase B, the global minimum point is equal to T,o{H) (see (|16|), (p!?!)), which is a /x-independent 
quantity. Hence, the particle density n = —dQ/dfi in the ground state of phase B is identically 
equal to zero, whereas in each phase Gk this quantity differs from zero. This conclusion follows 
from the definition of the thermodynamic potential fl given in the previous subsection , as well 
as from the fact that Q has no /x-dependence in the phase 

If fi increases, the curve of Fig. 5 moves up and to the right side of this figure. So, for 
some values of /i the function 0(E) in ( PT| ) will have three zeros (see Fig. 6): the one on the 
left-hand side in this figure, E^+i, is a local minimum of the function Vh^(E) , the one at the 
right-hand side, (S^ > S^+i), a global minimum of that function. If the chemical potential 
persists to grow, then at some critical value of /i the global minimum jumps from Sfc(yU,iJ) 
to Sfc+i(/x,if). At this moment we have a first-order phase transition from the massive phase 
Gk to the massive Gk+i one. In Fig. 4 the point of this phase transition lies on the curve 

Mtk+i, which is the boundary between regions Gk and Gk+i- Hence, all lines Mtk in Fig. 4 
{k = 0,1,2,...) are first-order phase-transition curves. Here we should also remark that all 

points of the line tQ^c{H) in this figure are described by Eq. (|23|) . Since the phase structure 
of the model I is so complicated, the dynamical quark mass S(/i, if), which is given by the 
global minimum of the potential Vh^^(S), also has a rather complicated /x, ii-dependence. For 
illustration, in Fig. 7 the schematic behaviour of S(/i, H) versus /i is presented at some fixed 
value of the external magnetic field H. 

From standard textbooks on statistical physics (see, e.g. [^]) we know that more than three 
curves of first-order phase transitions (lOPT) should not intersect at one point of a phase 
diagram; not more than three phases are thus allowed to coexist in nature. However, in Fig. 4 
one can see that infinitely many curves of lOPT cross at the one point M. Indeed, this visible 
contradiction with the above-mentioned statement is only fictitious, because at the point M 
we have a second-order phase transition |]T3. Hence, M does not belong to any of the critical 



curves of lOPT Mtk and is therefore not a point of phase coexistence. 



^ Since for all points of region B we have fi < So(-ff), it follows from (^) that Q — yf/^(I]o(iJ)) = (So(i^)), 
i.e. fl has indeed no /i-dependence. 
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4. Magnetic oscillations in the model I 



Now we want to show that there arise, from the presence of infinite sets of massless Aj. phases 
as well as of massive ones, magnetic oscillations (the so-called van Alphen-de Haas-type 
effect) of some physical parameters in the model I gauged by an external magnetic field. 



4.1 The case G < Gc 

Let the chemical potential be fixed, i.e. /i = const > /ic(H). Then on the plane {fi,H) (see 
Fig. 3) we have a line that crosses critical lines h, I2, ... at points Hi, H2, ... The particle density 
n and the magnetization m of any thermodynamic system are defined by the TDP in the 
following way: n = ~dQ/dfi, m = —dQ/dH. At /i = const these quantities are continuous 
functions of the external magnetic field only, i.e. n = n{H), m = m{H). We know that all the 
second derivatives of fl{fi, H) are discontinuous on every critical line /„. The functions n{H) 
and m{H), being continuous in the interval H G (0,cxd), therefore have first derivatives that 
are discontinuous on an infinite set of points Hi, Hk, ... Such a behaviour manifests itself as 
a phenomenon usually called oscillations. 

In condensed-matter physics [||, H it is a conventional rule to separate the expression for a 
physical quantity with oscillations into two parts: the first one is called monotonic and does not 
contain any oscillations, whereas the second part, which is of particular interest here, contains 
all the oscillations. Following this rule, we can write down, say, the TDP (^) of the NJL model 
I in the form 

0(/x, H) = n^onifi, H) + ^]„,,(/i, H). (29) 

In order to present the oscillating part V^osc{^J'■, H) as well as the monotonic one Q^onif^, H) in 
an analytical form, we shall use the technique elaborated in 0, where manifestly analytical 
expressions for these quantities were found in the case of a perfectly relativistic electron-positron 
gas. This technique can be used without any difficulties in our case, too. So, applying in ( pSf ) 
the Poisson summation formula 

00 00 °? 

^ a„$(n) =2^afc / ^{x) co?,{2'nkx)dx, (30) 

n=0 fc=0 

one can get for VLmon{l^^ H) and VLosc{l^^ H) the following expressions 

,4 ( ^TJ\2 



^mon 



^osc = -^y[^] \Q{Tiku)cosiT^ku + T^/A) + Pi-nku)cosi-nku -T^/A)], (32) 
' ^ \ 7ik J 



where u = ^ /{eH). (To find (31) and (R2[) it is sufficient to let the electron mass pass to zero 



in formula (19) of jg].) Functions P{x) and Q{x) in (|3TD and (|3^ are connected with Fresnel 
integrals G{x) and S{x) |]30| 



1 / X 

C{x) = - + i —[P{x) sinx + Q{x) cosx] 

2 y zTT 
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S{x) = \ — [P{x) cosx — Q{x) sinx]. 

They have, at x — oo, the following asymptotics pO|] : 

3 1 15 

P{x) = x~^ - -x'^ + Q{x) = --x"^ + —x~^ + ... 

4 ' ' 2 8 

Formula ( ^2]) presents, in a manifestly analytical form, the oscillating part of the TDP (^) for 
the NJL model at G < Gc- In the case under consideration, since the TDP is proportional to 
the pressure of the system, one can conclude that the pressure in the NJL model oscillates when 
i7 — > 0, too. It follows from ( |32D that the frequency of oscillations over the parameter (eH)^^ 
equals /i^/2. Then, starting from (|3^), one can easily find the corresponding expressions for 
the oscillating parts of n{H) and m{H). These quantities oscillate at H ^ with the same 
frequency and have a rather involved form, so we do not present them here. 

Finally, we should note that the character of magnetic oscillations in the NJL model at 
G < Gc resembles the magnetic oscillations in massless quantum electrodynamics ^. This 
circumstance is conditioned by the resemblance of the vacuum properties in the two models. 
Indeed, both in the NJL model and in QED, for fixed fi and varying values of H, there are 
infinitely many second-order phase transitions (see the appendix, where the vacuum structure 
of QED at /i, if 7^ is considered). 



4.2 The case Gc < G < (1.225...)Gc 

Here at /x > /iic(M) (see Fig. 4) the TDP of the NJL model as well as all thermodynamical 
parameters of the system oscillate with the frequency This can be shown in a way similar 

to what was done in the previous section. 

However, at M < /i < fiic{M), the character of magnetic oscillations is changed. Let us 
next prove this. First of all, in this case we have another expression for the TDP of the system 

nif^,H) = VH,inix,H)), (33) 

where Vh'^(E) is given in ( PP] ) and S(/i, H) is the non-trivial solution of the stationarity equation 
(pTf ). In each of the massive phases Gk the /x- and if-dependent function T,{fi,H) coincides 
with the corresponding fermionic mass Efc(/i, if) (see subsection 3.3). Using in ( PPD again the 
Poisson summation formula (|30|), one can easily select the oscillating part of the TDP (|33|) in 
a manifestly analytical form 

n,sc{l^,H) = 4^3/2 E (^J [Q(vrH cos(27rfccu + 7r/4) + 

+P{'Kkv) cos{2nkLJ - 7r/4)], (34) 

where u = /i^/(eif), u = {jj? — H))/{2eH). From (|3^ one can see that as a function of 
the variable {eH)~^, the TDP (^) oscillates with frequency (yU^ — S^(/i, H))/2 if this variable 
tends to infinity. Since f2(/i, H) is equal, up to a sign, to the pressure in the ground state of 
the system, also in the present case the pressure in the NJL model is an oscillating quantity. 
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Moreover, taking into account the stationarity equation (pTD, one can easily derive manifest 
expressions for oscillating parts of other thermodynamic quantities such as particle density 
n = —dQ/dfi and magnetization m = —dQ/dH from (|33|) and (0), 

ifi'^ - i7))3/2 ^ sm{2'Kkuj - 7r/4) 
m_ = -y^ + 



(/i^ - E2(/i, i7))3/2 - sm{2nkiu - n/A) 



S ' (35) 



It is clear from (|35|) that particle density and magnetization in the ground state of the NJL 
model oscillate with the same frequency as Q. 

Comparing magnetic oscillations in the NJL model I and in QED, we see three main dif- 
ferences. Let us remark that in QED the frequency of magnetic oscillations (over the variable 
{eH)~^) is equal to (/x^ — M^)/2, where M is the electron mass [^]. In the NJL model, in 
contrast to QED, the magnetic oscillation frequency is an if-dependent quantity. So, strictly 
speaking, in the NJL model magnetic oscillations are not periodic ones. Similar peculiarities 
of magnetic oscillations are observed in some ferromagnetic semiconductive materials such as 



HgCr2Se4 |^T|, where non-periodic magnetic oscillations over the variable {eH) were found 
to exist for electric conductivity^ as well as magnetization. This is the first distinction. 

A second difference is the character of the oscillations in the two models: in QED, magnetic 
oscillations are accompanied by second-order phase transitions (see Appendix), while in the 
NJL model they occur as a result of an infinite cascade of first-order phase transitions. 

Thirdly, we should remark that in the NJL model not only thermodynamic quantities os- 
cillate, but some dynamical parameters of the system do as well. This concerns, in particular, 
oscillations of the dynamical quark mass. In fact, by applying the Poisson summation formula 
(|30D in the stationarity equation (^1]) and searching for the solution Il{fi,H) of this equation 



in the form S(/i, H) = Smon + '^osc, one can easily find the following expressions for H 0: 

{fi^ - M^fl^ ^ sin(27rA;cu - 7r/4) 



y27r/icD3/2/'(M) A;3/2 Vcu3/2 

/i(/i2 - M2)3/2 



/i2-M2 
/i 



+ o[^,], (36) 



where u = (//^ — M^)/2{eH), M = So(/^) is the quark mass at H = 0, ft (see Fig. 2) and 



f{x) = F{x) + 2/iA/yu2 - x2 - 2x' In 



X 



(-F(S) is defined in (7)). Hence, in the framework of the NJL model I the quark mass S(/i, H), 
as well as other dynamical quantities composed from it, oscillate in the presence of an external 
magnetic field. 



^Magnetic oscillations of the electric conductivity, which is proportional to the particle density, are known 
as the Shubnikov-de Haas effect llsll. 
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5. Magnetic oscillations in the model II 



Next, let us consider magnetic oscillations in the more realistic NJL model II containing two 
kinds of quarks: u- and d-quarks with electric charges Ci and 62, respectively. The effective 
potential Vhfj. in the case under consideration is a trivial generalization of Vhh derived in the 
one-flavour case 

= -^-^ + ^2 (37) 

where S = Vcr^ + tt^, and ^./^^(S) is equal to Vh^{^) (pOf ), with e replaced by |ej|. 

Qualitatively, the phase structure of the NJL model II is the same as that of model I. So, 
at G < Gc we have an infinite set of massless phases (similar to the phase portrait of the 
model I in Fig. 3) reflecting the infinite set of Landau levels that is the basis for magnetic 
oscillations. Using the analytical methods of section 4, one can easily select in the present case 
the oscillating part of the thermodynamic potential: 

2 00 /|g.|//\^/^ 

= y-^j IQi'^^'^i) cosiirkui + tt/4) + P{Tiku,) cos{Ttku, - 7r/4)], (38) 

where z/j = fi^/{\ei\H). Hence, in the model II, in contrast to the model I and QED, we have 
a superposition of two oscillating modes. At growing values of the parameter {eH)~^, the 
frequency of oscillations in each of the modes is equal to e/i^/(2|ej|). 

At Gc < G there is a finite vicinity of Gc in which the phase portrait of the model II in the 
{fi, i7)-plane is similar to the one in the model I (see Fig. 4). So, in the case under consideration 
we have infinite sets of massless and massive phases as well. The cascade of massive phases 
is the foundation for non-periodic magnetic oscillations. Indeed, from (|37|) it follows that the 
oscillating part of TDP has the form 

^oscif^, H) = E E {^-^y IQi^ku,) cos{2nku, + 7r/4) + 

+P{'Kkui) cos{27rkuJi - 7r/4)], (39) 

where Ui = /x^/(|ej|if), Ui = (/i^ — S^(/i, H)) / {2\ei\H) , and S(/i, H) is the global minimum point 
of the effective potential (37). As in the previous model, E(/i, H) in the present case is an in- 
dependent function. This means that magnetic oscillations in the NJL model II are composed 
of two non-periodic harmonics, because each of them has, as a function of the variable {eH)~^, 
the if-dependent frequency e(/i^ — S^(/i, if))/(2|ei|). 

6. Summary and conclusions 

In the present paper we have studied the magnetic properties of a many-body system of cold 
and dense quark matter with four-fermion interactions. In particular, we have investigated the 
ground-state (vacuum) structure of two simple NJL models with one or two quark flavours, 
respectively, which are taken at non-zero chemical potential /z and magnetic field H. 
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As it turns out, in both types of models there exists a phase B (see Figs. 3 and 4) in 
which the quark mass is equal to So(if), i.e. it is a /^-independent quantity. Since this phase is 
achieved in the region /i < So(i/), the resulting particle density is expected to be zero, which 
is supported by our calculation. Clearly, this is in agreement with the physical interpretation 
of the chemical potential as the energy required to create one particle in the system. Indeed, 
the energy /i, which in the phase B is smaller than the quark mass, is not sufficient to create 
a particle, so that in the ground state of this phase the particle density must vanish. 

Most interestingly, we have shown that in NJL models there exist an infinite set of massless 
chirally invariant phases (phases Aq, Ai, ... in Figs. 3 and 4), which lead to periodic magnetic 
oscillations of some thermodynamic quantities of the system (so-called van Alphen-de Haas 
effect). In NJL models, this effect is observed at weak couplings {G < Gc, where Gc = 47r^/A^) 
or at sufficiently high values of the chemical potential and resembles magnetic oscillations in 
massless QED. 

Furthermore, for some finite interval of the coupling constant Gc < G < Gi, where in 
the framework of the model I Gi = (1.225...)Gc, the phase structure of NJL models I and II 
contains an infinite set of massive chirally non-invariant phases (phases Co,Ci,... in Fig. 4). 
This is the basis for non-periodic magnetic oscillations of some thermodynamic parameters, 
since the dynamical quark mass in each of the phases Gk is now itself an if-dependent quantity. 
(Notice that analogous non-periodic magnetic oscillations were recently found to exist in some 
condensed-matter materials too.) We should also remark that, at Gc < G < Gi, some 
dynamical parameters in NJL systems, such as quark masses, oscillate over [eH)~^ as well. 
This is an unknown fact in the standard condensed-matter theory of the van Alphen-de Haas 
effect. 

Moreover, our numerical analysis shows that for values of the coupling constant G G 
{Gi,G2), where G2 ~ 4:0Gc, we have in both models a finite number of massive phases Gk 
(the number of massless phases is infinite as before). At larger values of the coupling constant, 
i.e. at G > G2, there exist no massive phases in the phase structure of NJL models at all, 
except the trivial phase B. 

Finally, it is interesting to note that for fixed magnetic field H the dynamical quark mass 
T,{fi,H) of NJL models discontinuously jumps as a function of /i, at points fiQ, fii, fi2, ... (see 
Fig. 7), thus refiecting the structure of the underlying phase portrait shown in Fig. 4. 

In conclusion, we have shown that NJL models at /x, if 7^ exhibit an interesting phase 
structure and a set of oscillating quantities, which are richer than in the corresponding QED 
case. Note also that the above approach can be further employed for a combined study of quark 
and diquark condensates of NJL models taken at fi,H 0. Work in this direction is under 
way. 

One of the authors (D.E.) gratefully acknowledges the kind support and warm hospitality 
of the colleagues of the Theory Division at CERN. This work was supported in part by the 
Russian Fund for Fundamental Research, project 98-02-16690, and by DFG-project 436 RUS 
113/477. 
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Appendix 



The QED vacuum structure in the presence of /i and H. 

In the framework of QED and in a one-loop approximation, the effective Lagrangian in the 
presence of an external homogeneous magnetic field has the following form [|, ^: 



Le//(/i, H) = Li{H) - nQEoilJ', H), 



(Al) 



where Li{H) is the Lagrangian at /i = 0, 7^ 0. Since it does not infiuence the phase structure 
of QED, we do not present its explicit form here, but refer to |^. The second part of 
Leff{fi,H) is exactly the thermodynamic potential of a perfect electron-positron gas 



{fi, H) = -— ^ aJifi - e„){/iyy 



QED[ 



47r2 



In 



n=0 



^2 - el + fi)/er. 



{A.2) 



In (A. 2) we have introduced the notation: e„ = V + 2eHn, where M is the electron mass 
at = 0, /i = 0, and a„ = 2 — 5on- 

Using the results and methods of subsection 3.2 one can now easily show that, on each line 
/„ = {{jji.H) : jjL = \/ M"^ + 2eHn} of the (/i, if)-plane, all second derivatives of Leff{fx,H) 
(A.l) are discontinuous. So, the lines In [n = 0, 1, ...) in Fig. 8 are the curves of second-order 
phase transitions. They divide the (/x, if)-plane into an infinite set of regions C„ (n = 0, 1, ...) 
corresponding to different massive phases of QED (associated to Landau levels) with the same 
electron mass M. Nevertheless, each phase C„ is characterized by such physical quantities 
as particle density n{fi, H) and magnetization m{fj,,H). On each phase boundary /„ these 
quantities are continuous functions. However, their first derivatives are discontinuous functions 
on each line In, so the derivative jump of n or m is the signal of a second-order phase transition. 

As in the previously considered case (see subsection 3.3 and Fig. 4) the point M in Fig. 8, 
where all lines /„ intersect, is a special point differing from other points of the lines /„. Indeed, 
at H = we have HI 



^QEoif^, 0) 



e{fi - M) 



M 



M2 



(/i 



X] 



'(/i + 2x). 



(A3) 



It follows from (A. 3) that only the third derivative d^VLQEoil^,^) / {d^Y is a discontinuous 
function at the point /i = M. At the same time, in other points of critical curves /„, already 
the second derivative d'^flQEoifJ', H) / (dfi)"^ is discontinuous. 

Here we have considered only the relativistic case, but one can easily show that the non- 
relativistic electron gas at /x, i7 7^ has an infinite cascade of massive phases, too. So, at 
the basis of the van Alphen-de Haas and Shubnikov-de Haas effects lies an infinite set of 
second-order phase transitions. 
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Figure captions 



Fig. 1 Phase portrait of the NJL model at non-zero // and for arbitrary values of the fermion 
mass M . Phases B and C are massive and non-symmetric, A is a chirally symmetric phase. 
Here ^2c{M) = M, fiic{M) = ^J^M^ \n{l + A"^ / M^) , HadM) is the solution of the equation 
V^{0) = V^{M), = A/(2.21...), Mic is the solution of the equation fij^iMi^) = AV(4e). 
In phase B the particle density in the ground state is equal to zero, whereas in phase C it 
is non zero. Solid and dashed lines represent critical curves of second- and first-order phase 
transitions, respectively; a and (3 denote tricritical points. 

Fig. 2 The behaviour of the dynamical quark mass So(/x) as a function of /j, for the case 
M < Mic and H = within the framework of the model 1. 

Fig. 3 Phase portrait of the gauged model I at G < Gc- Solid hues Ik are given by 
k = {if^,H) : ^ = V2eHk}. They are critical curves of second-order phase transitions. The 
dashed line of first-order phase transitions is defined by Eq. (23). 

Fig. 4 Phase portrait of the gauged model I at Gc < G < (1. 225. ..)Gc. Here M is the quark 

mass Bit fj, — 0, H — 0, and the quantity /Xic(M) is presented in Fig. 1. The dashed line tofidH) 
is defined by Eq. (23). In this case one has infinite sets of symmetric massless phases Aq, ^i, ... 
as well as massive phases Go, Gi, ... with DCSB. In addition there exists another massive phase 
B. 

Fig. 5 Typical behaviour of 0(S) (21) for some points (/x, H) E C^- Here (7„ = 
Va*' - 2eHn. 

Fig. 6 For increasing there arise three zeroes of the function 0(S) (21) defining two local 
minima and Efc_|_i(Syfc_|_i < S/j) of the effective potential, respectively. The global minimum 
of ^///i(S) lies in one of them and passes by a jump from one local minimum to another one 
depending on the values of /x. 

Fig. 7 Schematic representation of the dynamical quark mass S(/i, H) as function of ji for 
fixed magnetic field qualitatively refiecting the structure of the phase portrait in Fig. 4. Here 
the magnetic field is fixed in the interval {Hs^,Ht^), where S2,^2 are tricritical points (see Fig. 
4). The ijLk,pi'c{H) (A; = 0, 1, 2) are the values of the chemical potential at which the line H — 

const crosses in Fig. 4 the critical curves Mtk and S2t2, respectively. 

Fig. 8 Phase portrait of QED ai ji^H ^ 0. All the fines Z„ are the curves of second-order 
phase transitions. 

Table caption 
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Table 1. Values of the external magnetic field corresponding to tricritical points Iq and Sq 
(see Fig. 4) for different ratios of coupling constants G/Gc- 
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G/Gc 


1.01 


1.1 


1.15 


1.2 


TT / A 9 


0.0129... 


0.08119... 


0.10769... 


0.12987... 




0.00614... 


0.05639... 


0.08088... 


0.10338... 
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